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PROBLEM 


Apply statistical decision theory to the design of automatic radar detection 
systems. The specific phase of the problem reported here is the theoretical investi- 
gation of a system based on a two-stage binomial test which uses low resolution 
data in the first stage and high resolution data in the infrequently employed 
second stage. 


RESULTS 


1. Detection probabilities for binomial two-stage detection systems with 
coarse-fine range resolution have been calculated; they indicate substantial 
saving over conventional systems. 


2. The use of sequential rather than fixed-sample-size testing in the first 
stage has been found, in the special cases analyzed, to improve the detection 
probability only slightly. 


3. In certain clutter-free situations, doubling the number of first-stage range 
bins used to cover the detection zone (by halving the first-stage pulse length and 
doubling the instantaneous power) has been shown to cause a loss in detection 
probability equivalent to roughly 0.3 to 0.5 dB/pulse in signal-to-noise ratio 


(assuming that the dwell-time, false-alarm rate, and second-stage range resolution 
are held constant). 


4. The optimum length of the second stage has been shown to vary consider- 
ably with other design parameters and quantities. 


RECOMMENDATIONS 


1. Extend the analysis of two-stage systems with fixed-sample stages to 
cases where the radar cross section is slowly fluctuating. 


2. Calculate the cumulative detection probability (the probability that an 
approaching target is detected before it reaches a given range) for cases of 
practical interest. 


3. Consider two-stage systems for applications in which the power capability 
of the radar dictates that the average number of pulses per beam position be large. 


4. Consider two-stage systems having both doppler and range resolution. 


5. Investigate the effect of clutter and interference on the performance of 
the two-stage system. 


6. Implement a two-stage system as soon as a suitable antenna is available; 
perform two-stage and conventional tests with live targets to verify performance 
predictions. 
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INTRODUCTION 


A theoretical analysis was made of the performance of some binomial two- 
stage radar detectors which use low range resolution in the first stage and high 
range resolution in the second stage. In the first stage of the detection procedure, 
one or more pulses are transmitted while the antenna beam is in a given position, 
and a tentative decision is made for each of the wide range bins as to whether or 
not a target is present. If no alarm occurs in the first stage, the antenna beam is 
moved to the next position; otherwise, the beam is held fixed and pulses of the 
second-stage type are transmitted. In the latter case, the regions corresponding 
to first-stage alarms are then re-examined by collecting and processing the high- 
range-resolution data, after which the antenna beam is moved to the next position 
for a new first stage. 

The two-stage detectors reported on here belong to a class of multiple-stage 
radar detectors requiring beam agility and versatile transmission and reception 
capabilities. After a brief description of this class of multiple-stage detectors and 
a mention of particular versions considered by others, the binomial high-resolution 
two-stage detectors are investigated and their performances compared with that of 
a conventional detector. 


HISTORICAL BACKGROUND 


The large body of literature on sequential analysis reveals a wide variety 
of hypothesis-testing situations in which a statistical test of the sequential type 
will require fewer observations, on the average, than a fixed-sample-size test with 
the same probabilities of error. With the development of electronic scan antennas, 
the application of sequential testing procedures to radar detection has become 
practical; in many detection situations it is possible to design sequential pro- 
cedures which promise appreciable energy savings over procedures limited by the 
constant rotation rate of a mechanically scanning antenna. 

In the study of sequential radar detection the usual assumption has been 
that successive observations (detected returns from a particular range gate) are 
obtained under nearly identical circumstances (identical transmitted waveforms, 
same receiver characteristics, etc.) and the testing procedures considered have 
typically consisted of operating on the acquired sequence of measurements of a 
matched-filter receiver output with the same rule of procedure after each trans- 
mission. In particular, earlier work was generally with Wald’s* sequential prob- 
ability-ratio test of two statistical hypotheses, the optimum nature of which can 
be realized in cases with a single resolution cell (or a single decision) per 
antenna-beam position. Later applications to radar involved modifications of the 
classical Wald test and various schemes for using sequential tests in multiple- 
resolution element situations. 


‘See list of references at end of report. 


Recently some attention has been given to a type of sequential detection 
in which the statistical testing procedure and the radar transmission and reception 
characteristics are programmed to change between pulses or groups of pulses. 

The general form of this type of muitiple-stage detection is outlined in figure 1. 
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Figure 1. Multiple-stage detection with a maximum of m stages. 


An electronically scanned antenna points, once each on every scan, in a number 
of different directions. For any given stay of the beam, stage-(7+1) testing occurs 
only if stage-i testing indicates that a signal might be present. The antenna beam 
stays in a given position until a final decision is made as to the presence of any 
targets in that sector. Note that in this type of multiple-stage test, final signal 
decisions can be made only in the final stage, while a final decision that no signal 
is in any resolution element can be made in any stage. Some of the many possible 
degrees of freedom in the design of each stage are the following. 


1. Range resolution, number and placement of range bins 


2. Doppler resolution, number and frequency coverage of filters 


3. Energy per transmission 


4. General class of statistical testing procedure (e.g., fixed-sample-size 
tests, Wald sequential tests,’ etc.) 


5. Number or average number of transmissions 

6. Choice of not quantizing data or quantizing data into any number of levels 
7. Decision statistic 

8. Decision thresholds and, if used, quantization thresholds 


A two-stage detection system having the same radar transmission 
characteristics in both stages was investigated by Helstrom.? In the system he 
describes, square-law detected returns from a number of pulses would be processed 
by superposing them and comparing the accumulation at the end of each stage with 
the threshold of that stage to see if it exceeds it for any delay, that is, for any 
range. The second-stage decision would be based on both first-stage and second- 
stage returns. For convenience in analyzing the process, he supposes that each 
return is sampled regularly in range, with the number of samples, or increments, 
equal to twice the product of the bandwidth and the sampling period. Letting the 
test statistic for each of the range increments be the sum of sampled returns over 
successive pulses is then analogous to integrating video returns. Two of 
Helstrom’s conclusions as to the performance of this system were that its advan- 
tage over the ordinary detection system (a fixed-sample single-stage system) 
decreases as the SNR increases and that at best it achieved about twice the 
average scanning rate attained by an ordinary detection system. His investigation 
was not extended to cases with three or more stages since it appeared that little 
further improvement could be expected and the computations would be much 
more difficult 

H. M. Finn analyzed an ‘‘Energy Variant Sequential Detector’? (EVSD)2° 
and, in a later report,° a ‘‘Resolution-Variant EVSD.’’ For rapidly fluctuating 
targets (Swerling’s’ fluctuating target types 2 and 4) the decision statistic at 
each stage is the sum of square-law detected outputs, generally for three pulses 
in the first stage and six in the second, while for slowly fluctuating (types 1, 3) 
and nonfluctuating targets a linearly detected matched-filter output is used, 
generally with two or three stages of one pulse each. In some of the proposed 
versions of an EVSD the sum of both first-stage and second-stage outputs is used 
in the second stage. References 4 and 6 include analytical results on the optimi- 
zation and performance of some of the EVSD’s described, a verification of pre- 
dicted performance by simulation on a digital computer, and descriptions of the 
general implementation of EVSD’s and RV-EVSD’s. One of his conclusions was 
that in many applications the increased efficiency which would result from 
employing more than two stages would not compensate for the inconveniences 
in implementation. 

A more thorough analysis of an EVSD having single-pulse stages is made 
by Brennan and Hill.*? They present optimum energy levels and cumulative 
detection probabilities for various false-alarm probabilities and numbers of range 


bins, and for each of the five standard target types. Their results indicate that 
the energy-variant two-stage radar has a power saving of several dB over a con- 
ventional, uniformly scanning radar. 

In a brief note in an IEEE Proceedings, R. F..Baum gives a comparison 
of some two-stage and three-stage tests and of a sequential probability-ratio 
test with a conventional detection mode.*° 


BINOMIAL HIGH-RESOLUTION TWO-STAGE DETECTION 


Glossary of Symbols 


Ma number of stage-A range bins per beam position 
Mp number of stage-B range bins in a stage-A range bin 
m M4 Mg, number of range bins searched by the single-stage system 


used as a comparison 


R range to target 
[Gass range to outer circumference of surveillance zone 
y ; 
r /R yay, normalized range 
S mean power signal-to-noise ratio (SNR) per pulse, expressed 


in decibels 


Sa SNR in dB/pulse corresponding to a ‘‘design’’ target at range R,,,., 
S(r) SNR in dB/pulse as a function of target range, with S(1) = S, 
S4 SNR in dB/pulse for stage A when it differs from that for stage B, S 


then denoting the latter 
FAR false-alarm rate (average number per hour) 
Pp, Pp(S) _ single-scan detection probability 


PREF pulse repetition frequency — number of pulses used per second by 
the two-stage mode only, if time-shared with other modes 


a, 8 parameters of Rice and Rayleigh distributions, respectively 


N, Ng, Np number of pulses in a stage 


N 


N(S) 


E(n), E,(n) 


oH YN Che 


B, Ba, Bp 


9, Wa, IB 


Do 
Pa, PB 


p(S), pa(S), 
PB (S) 


K, Ky, Kp 


C 


average number of pulses per beam position when noise alone is 
present, for a two-stage detector 


average number of pulses per beam position occupied by a target, 
for a two-stage detector 


average number of pulses in a sequential stage-A, E,(n) in the 
noise-alone case 


false-alarm probability per test — same for every range bin in that 
stage; a, is stage-A false-alarm probability when M4 = 1 


miss probability per test (therefore per stage) 


quantization level against which sample voltage x is compared; 
level may vary with range 


prob {x > q} for the sample voltage x corresponding to any particular 
stage and range bin 


value of p when noise alone is present 
specific cases of p, 


value of p for particular values of the SNR 


decision threshold of a fixed-sample-size test 


parameter specifying a sequential test 


General Description 


The detection scheme to be investigated here is a binomial two-stage 
system featuring coarse resolution in stage A and fine resolution in stage B. The 
properties of the stages are listed below. Binomial quantization is performed by 
comparing each of the sample voltages with a quantization level chosen for that 
particular stage and range bin. 


STAGE A 


1. Long pulses, generally few range bins. The long pulses result in radially 
wide range bins; thus the surveillance zone is divided in range into a relatively 
small number M , of bins per beam position. 


2. One stage-A decision per stage-A bin per beam position. There are M4 
stage-A tests per beam position, each having the two possible decisions, alarm 


and no alarm. Computations will be for both binomial sequential testing and 
binomial fixed-sample-size testing. 


3. Very few pulses. In a given beam position, a small number N, of pulses 
are transmitted for stage-A fixed-sample-size testing, or a variable number for 
sequential testing. E, (n) denotes the mean number of pulses in a sequential 
stage-A when only noise is present. 


STAGE B 


1. Stage B used infrequenily. In a given beam position, a stage B is used 
only if one or more alarms occur in stage A. 


2. Short pulses or pulse compression, many range bins. The My wide bins 
of stage A are each subdivided into Mp radially narrow bins. (The equivalent 
single-stage system must therefore examine M, Mp bins.) The term pulse will 
refer to the entire waveform transmitted each pulse repetition period. 


3. Not all bins examined. Stage-B tests are performed for those stage-B 
bins contained in a stage-A bin which alarmed.* Since the probability of a false 
alarm in stage A is small, the probability of multiple false alarms in stage A is 
very small; thus the number of stage-B bins examined when a stage B is used is 
usually Mp in noise-only cases. 


4. Same pulse repetition frequency as in stage A. 


5. Many pulses. Stage-B testing is the performance of a number of binomial 
fixed-sample-size tests, one for each bin being examined. The sample size Np of 
these tests is considerably larger than N y or E,(n). 


6. Same SNR as for stage A in results shown. However, the equations given 
throughout are also valid when the SNR (signal-to-noise ratio per pulse) is dif- 
ferent in stage B. 


Distribution of Signal and Noise 


The graphical results presented assume pulse-to-pulse independence and 
are for the Rayleigh and Rice cases, which correspond respectively to Swerling’s 
rapidly fluctuating target type 2 and nonfluctuating target type. Computations for 
the rapidly fluctuating target of type 4 were made but for brevity are not included, 
since in every situation the type 4 curve fell almost entirely between the curves 
for type 2 and the nonfluctuating type. 


*It would be better to examine, for each stage-A bin that alarms, the stage-B bins contained 
in a somewhat larger range interval (covering the alarming stage-A bin). An analysis of this 
more complicated case is unnecessary, since the higher threshold(s) required to provide the 
same false-alarm rate generally result in only a very slight loss in detection probability. 


For target types with only scan-to-scan independence (slowly fluctuating 
types | and 3 for example) it would be necessary to caiculate the single-scan 
detection probability 


PL )= feo (5:5) dS 


where g(S;S) is the density function of the SNR with mean S and Pp (S) is the 
single-scan detection probability for a given value of the SNR. 

The equations below are stated in the form for linear detection, but the 
results based on these may be used for square-law detection as well. 

When the output voltage x of an envelope detector has the Rayleigh 
probability density function 


f(x) = — exp (-x?/20?) (x > 0, 20? = E(x?) 


oO 


where o*= 1 when only noise is present and o” > 1 when a signal is present, 
then the probability that x exceeds a quantization level q is 


PD = Po = exp(—q’/2) 
in the noise-only case and is 

p = p(S) = exp(-q’/20’) 
in the signal case, where 

S= 10 log,, (207 ~1) 


is the mean power SNR in dB/pulse. 

For a nonfluctuating target the normalized output voltage x (in units of 
rms i-f voltage) after linear envelope detection is assumed to have the well-known 
Rice density function 


f(x) = x exp ez = Jota) 


where a is the ratio of peak i-f signal voltage to rms i-f noise voltage and I,(+) is 
the modified Bessel function of first kind, zero order. The probability that x 
exceeds gq is 


P = Po = exp(-q’/2) 


in the noise-only case and is (ref. 7, page 159) 


= p@= on ESS (q/a) 1 (aq) 
jel 
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in the signal case, where 
S= 10 log,, a’/2 


is the mean power SNR in dB/pulse. 

Let S, be the SNR corresponding to a ‘‘design’’ target (a target with 
specified radar cross section) at range R_... The fourth-power relation between 
range and received echo-power gives as the SNR of this design target at range r 


S= S(r) = Sz 0 log,, r 


where r is the normalized range R/R..,. 


Binomial Fixed-Sample-Size Testing 


In a stage of the fixed-sample type, a binomial fixed-sample-size test is 
performed for each range bin. An individual test is performed by comparing, for 
each of N successive independent pulses, the sample voltage (the amplitude of 
the envelope at the time delay corresponding to the given range bin) with the 
quantization level qg, and deciding ‘‘alarm’’ if at least K of the N voltages exceed 
q and ‘‘noise only’’ otherwise. If p is the probability for a particular range bin 
that the sample voltage for any given pulse will exceed the threshold q, the 
probability for that range bin that K out of N independent samples exceed q is the 
cumulative binomial probability 


N 


D(N,K,p) = » (%) pi(1—p)N-i 
i=K 


(For convenience the subscripts designating stage are omitted here.) For fixed N 
and a each possible value of the threshold K determines a value of p,, from 

a= D(N,K,p,), and therefore determines a value of g. Then each K, from its q, 
also determines a value of p(S) for any given value of S. The optimum value of K 
for the given value of S is that which minimizes the miss probability B(K,S) = 

1 —D(N,K,p(S)). Whenever convenient, 8(K,S) will be written B(S) or 8, and 

since the optimum values of the thresholds are used for each S in all calculations, 
B(S) and 8 will be the minimum B(K,S) for that particular N, a, and S. Ideally, for 
each range-bin test the threshold K should be chosen to optimize detection for 
the SNR of the design target at that range, i.e., for S(r), where arbitrarily r could 
be the center of the bin. Generally, however, the use of a different quantization 
level q for each bin is impractical in implementation. It is fortunate, therefore, 
that using the same value of K for all bins results in miss probabilities which do 
not differ appreciably from those in the ideal case. In figure 2, curves of 8(K,S) 
are given for a typical stage-B test (N= 12; a=6 x 107°; Rice distribution); note 
that K=7 is optimum over effectively the entire range of 8(S), as S varies, and 
that the curves are relatively flat over several values of K. 
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Figure 2. These curves for a binomial fixed-sample-size test typical of those 
employed in stage B indicate that using the same threshold for all range bins 
would result in little or no loss in optimality. 


False-Alarm Rate (FAR) 


For a fair comparison of the detection capabilities of various testing 
procedures, all should have the same FAR. The FAR of the fixed-sample single- 
stage system is proportional to the number m of range bins per beam position and 
to the false-alarm probability a of an individual test and is inversely proportional 
to the number N of pulses per beam position. 


FAR = (3600 sec/hr) x PRF x m a/N 


The FAR of the two-stage system (sequential or fixed stages) has no simple 
exact formulation but is closely approximated by 


FAR = 3600 x PRF x M, Mp a,ap/N 
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where «,ap is the per-bin false-alarm probability for all M, Mp range bins and 
N is the average number of pulses per beam position when only noise is present. 
When comparisons are made where the number of range bins and the PRF are 
fixed, the FAR is made the same for all cases by holding «,ap constant if N is 
fixed or a4ap/N constant if N varies. 


Sample Sizes 


When applying two-stage detection to a particular detection situation, it 
probably will not be immediately evident what values of N, N, (or E,(n)), Np 
and a, are best to use. A relation among these quantities is 


——_ MA 

N= Na+ 1-(1-ay) Np 
for fixed-sample stages, and 

N= E,(n) + a,Np 


when the first stage is sequential and involves a single range bin. (Sequential 
tests are discussed later; calculations were not made for the multiple-resolution- 
element type.) 

The nature of binomial tests makes it necessary to use numerical rather 
than analytical methods to find the optimum values of the unspecified quantities, 
but the amount of computation time required on a high-speed computer is moderate. 

When values of N, My, Nq (or Ep(n)), S-S,,* S, and ayag (or «,ap/N) 
are specified, the optimum Np is the value yielding the maximum single-scan 
detection probability 


P,(S) = [1-2.4(8.)| i-en(S) 


For each value of Np considered, the values of «, and ag must be determined 
before calculating Pp(S); the value of «4 is determined from the equation for N 
(an additional relation between a, and E,(n) must be satisfied in the sequential 
case), and ap is found from fixed a,ap. It is important to keep in mind that the 
optimum value of Np varies with S and therefore varies with range for a target of 
given cross section. For simplicity it is assumed that the same number !Vp of 
pulses are transmitted no matter which stage-A bin alarms. 

The optimization of N,, given N, and the optimization of N are 
discussed later. 

Another sample size of interest is the mean number of pulses when a 
signal is present. When M,= 1, this quantity is N(S) = N4 + [1-f(S)] Np for 
the fixed-sample case and is N(S) = E(n) + [1-p(S)} Np for the sequential case, 
where E(n) is a function of S. 


*S, and S are the per-pulse SNR’s in stage A and stage B, respectively; their difference is 
assumed constant with range and cross section, and is zero for all results presented. 


Resolution 


Large range bins for stage A can be produced by transmitting pulses of 
long duration. When only one range bin is used in stage A, the minimum range 
of the bin must be at least half of R,.,, because of propagation time. The two- 
stage mode would then be an acquisition mode, and a single-stage high-resolution 
mode could be included for close targets; in fact, for any application in which each 
stage-A pulse is very long, detection for the resulting blind range must be provided 
for in some other way. If each stage-A pulse is a very long uncoded waveform, the 
width of the signal frequency spectrum might be small compared with the doppler 
frequency shift. Then, to cover the range of expected doppler frequency, it would 
be necessary to use in stage A either a bank of doppler filters or a receiver band- 
width considerably larger than the widthof the signal frequency spectrum. The 
analysis of a two-stage system with doppler filters used in this manner is more 
complicated and will not be considered here. If the doppler shift is provided for 
by making the receiver bandwidth large or if sampling long pulses results in a 
loss in SNR, then the decrease in the stage-A SNR must be taken into account. 
As the comparisons presented in this paper are in terms of a per-pulse SNR 
common to both stages, power saving can be measured only when assuming that 
both the transmitted energy and the resulting SNR are the same for stage-A 
pulses as for stage-B pulses. Otherwise, one must compute for each case the 
average amount of transmitted energy required per beam position. 

Even though the resolution cells will be referred to as range bins, the 
results hold for any type of resolution such that (1) each of the My cells is sub- 
divided (in range and/or doppler frequency or other) into Mp cells, and (2) a 
stage-B cell is tested if and only if the one stage-A cell containing it alarms. 

Accepting the fact that the resolution-varying property of two-stage 
detection contributes to its efficiency, the logical conclusion is that the range 
of surveillance (in a clutter-free environment) should be covered with as few 
stage-A range bins as possible under the limitations on minimum range and 
doppler shift mentioned above. This is supported by the Pp-versus-S curves 
given in figure 3 for fixed-sample stages, a Rice distribution, and N= 1.2. 

Since the optimum value of Ng decreases with increase in M,, each curve was 
computed using the Np which was optimum for Pp about 0.5; these values of Ng 
were 14, 10, 7, 5, and 4, respectively, for M,= 1, 4, 16, 64, and 256. Note that 
quadrupling the number of stage-A range bins causes an increase of roughly 0.6 
to 1.0 dB/pulse in the SNR required for a single-scan detection probability of 
0.5. For a fair comparison, the M, bins would cover the same search zone in 
each of the cases; then the product M, Mp must be the same for each case in 
order that the FAR and final resolution be the same. (FAR = 3 PRF x M4 Mpg x 10~ 
The correspondence of range bins (bin centers) to values of S is determined by 
whatever value is assigned to S,. The Pp-versus-S curves are valid for any Sy 
(in the range of the graph) as long as the threshold assigned to each range-bin 
test is the optimum one for the value of S corresponding to that range (for that Sy) 
and the target is indeed of the design cross section. For any other situation, 
where S is the actual rather than design value of the SNR, the curves give only 
upper bounds to Py, generally close bounds, though, since (as was suggested by 
fig. 2) the thresholds would still be at least nearly optimum. 
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Figure 3. Single-scan detection probability vs required SNR, where the final 
resolution remains constant but the coarseness of the stage-A bins varies. 
Rice distribution; N = 1.2, optimum Np for eachM, for Pp = 0.5; FAR = 

3 PRF x My Mp x 10~, where M, Mp may vary. 


The case of a single range bin in stage A will be treated in detail since 
(1) the detection zone to be searched by a two-stage system will in some applica- 
tions be sufficiently narrow radially that the pulse length necessary to bridge it 
with a single range bin is not inefficiently long, and for these applications M, = 1 
is optimum; (2) from the analysis of the M, = 1 case, many general conclusions 
can be drawn concerning the nature of two-stage systems having any M,; and 
(3) 1t is the simplest case to analyze. 


N and the Average Scanning Rate 


The average scanning rate in revolutions per minute for a two-stage 
system is 60 x PRF / (no. of beam positions) xN when no target is present and 
is approximately that when the target density is low. Suppose first that the PRF 


is predetermined by equipment restraints or other considerations. The following 
properties of the two-stage system are pertinent to the problem of selecting a 
suitable value of N. 


1. The single-scan detection probability Pp increases with N. Figure 4 
illustrates this and also indicates how the optimum value of Nj, changes with N. 
(Fixed-sample stages; Rice distribution; S= 5 dB/pulse; M, = 1; asap = Nx 10° ae 
K,4, Ng, Kg optimum; FAR = 3.6 x Mg x PRF x 1075, where Mp may vary.) 


2. The average number of scans across a target staying in the surveillance 
zone a relatively long time is approximately inversely proportional to N if the 
number of beam positions is large and the target density is low. 


3. The cumulative probability of detection (probability of at least one detec- 
tion before a given range) tends to increase slightly with N but in an irregular 
manner. Property 3 is the resultant of properties 1 and 2. It has not been verified, 
but is supported by calculations of the quantity 1-(1 -Pp yconstant/N where Pp is 
a function of N and of the optimum test parameters for that N. This quantity is an 
indicative measure of the cumulative detection probability; it is based on the 
fact that the probability of detection is equal to 1+(1- Pa for a target scanned 
L times while at the same range, with L approximately inversely proportional 
to N by property 2. Curves of 1-(1-P maa versus N are given in figure 5 for 
the values of Pp plotted in figure 4. 


4. The length of an optimized two-stage test in the case of a detection (or 
at least a stage-A alarm) tends to increase with N. Figure 6 illustrates this 
with a plot of Ny +Np versus N, where Nz is optimized for each Ny and N 
Although N , + Np is an upper bound to N(S), it varies with N in a similar 
manner and is approximately equal to N(S) for signals somewhat stronger than 
that for which the Np’s are optimized here. 

Properties 3 and 4 suggest that a small N is at least nearly as efficient 
as a large N, so that choosing the most suitable average scanning rate is largely 
a matter of determining (and this involves examining the corresponding single- 
scan detection probability curve) the most desirable value of the number of scans 


across a target (of particular velocity) while it travels through the surveillance zone. 


Consider now the case where the PRF is adjustable in that it may be set at 
one of a wide range of values. The optimum choice of the combination of PRF, 
average scanning rate, and N then depends on the particular application. For a 
situation in which the radar is peak-power limited and the cost of the power can 
be disregarded, the PRF might be best set at the highest rate allowing the de- 
sired maximum unambiguous range. If, on the other hand, the major limitation is 
that of average power, then the most efficient procedure would probably be to 
maximize the transmitted energy per pulse by using a small N and, for that N, the 
lowest PRF that permits the antenna to scan across the target a reasonable 
number of times. 

It is fair to conclude that a small N is suitable in many applications. Also, 
note in figure 4 that for N as large as 2.6 the Ny, = 1 case is still better (in the 


particular situation considered and in other typical situations) than the N A = 2 case. 


The subsequent analysis is restricted to cases with Ny = 1 and E,(n) < 2 for these 
two reasons and because computationally they are the most tractable. 
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Figure 4. Single-scan detection probability vs mean number of pulses (noise- 
only case) per beam position. Rice distribution; S= 5 dB/pulse; My, = 1; 
K,, Ng, Kp optimum for each point; FAR = 3.6 Mp x PRF x 10™°, where Mp 


may vary. 
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Figure 5. An approximation to the probability of detecting a constant-range 
target when the mean number of scans across the target is 4/N; an indication 
of how the cumulative detection probability would vary with N. The values of 
Pp used are those plotted in figure 4. 
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Na + optimum Ng 
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NOTE: Na=2 case 
dips when optimum Ka 


changes from 2 to |. 


Figure 6. Number of pulses in an optimized two-stage test when a target is 
detected vs the mean number when noise only is present, for the cases plotted 


in figures 4 and 5. 


ANALYSIS OF SPECIAL CASES 


General Description 


Sequential Cases. STAGE A: Binomial Wald’ sequential test, single- 
range-bin type (M, = 1), with E,(n) less than about two. STAGE B: Mp range 
bins (radially narrow). Returns from Np pulses used in Mg binomial fixed-sample- 
size tests. 


Fixed-sample Cases. STAGE A: One pulse, one range bin (N,=1, My = 1). 


STAGE B: Mp range bins. Returns from Ng pulses used in Mp binomial fixed- 
sample-size tests. 

The detection capabilities of these two-stage detectors will be compared 
with that of a conventional single-stage detector, one that employs a binomial 
fixed-sample-size test for each of m= Mp range bins. 


A Sequential Test For Stage A 


When there is only one resolution element to be tested and E,(n) is to be 
less than about two, the ‘‘best’’ binomial sequential test is a very simple form 
of Wald’s test, one which may be performed in the following manner. A random 
variable y; is assigned the value 1 if the sample voltage x for the ith pulse exceeds 
q, and the value (-1) otherwise. After each pulse, say the nth, one of three 
possible actions is taken. (C is an integer larger than two and is specified in 
accordance with desired test characteristics.) 

If0 <1+y,+---+y, <C, testing is continued with another pulse. 

If C <1+y,+++-+y,, testing is terminated with an alarm. 

If l+y,+-+-+y, <0, testing is terminated without an alarm. 

A test of this form is ‘‘best’’ in the following sense. 


1. It is a probability-ratio test for (and only for) any p, and p(S) such that 
p(S) = 1-p,. (The probability that y; equals 1 is p, when noise only is present 
and p(S) when the SNR has value S.) Since the SNR varies widely with range and 
cross section, the choice of the value S for which p(S) = 1—p, as the value for 
which to optimize the test is often as fair as any. 


2. For any binomial sequential test having an E,(n) < 2, there is some C 
and q such that a test of the above form has the same E,(n), approximately the 
same a, and a B(S) which is not appreciably larger for any S (assuming any 
particular functional form of the noise and signal-plus-noise distributions and 
assuming q adjustable for each test). This was apparent from an examination, 
for a wide variety of binomial sequential tests, of the trade-offs in detection 
characteristics that occur from test to test; the nature of binomial tests does 
not permit direct comparisons. 


3. A test of this form is much simpler to analyze than any other binomial 
sequential test. Note that the sum 1+¥y +--++y,, will equal one of the bounds 
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at termination. Since Wald’s approximation equations yield true values for tests 
in which there can be no excess over either bound, exact results can be obtained 
without resorting to unwieldy iterative computational methods. The necessary 
equations can be put in terms of quotients of polynomials;* those for the C = 6 
case, for example, are 


5 


2 4 
E(n|p) = 1—2p + 3p" +3 , Ptalarm|p} = es 
1-4p +7 p?-6p?+3p* 1-4p +7p*—6p? + 3p* 


E,(n) = E,(n|p,) is assigned a value, and the computer solves for p, =p, and 
a, = P{alarm|p,} and, given a value of S, for p(S), B4(S) = 1—Ptalarm|p(S)}, and 
and E(n|p(S)). 


Detector Optimization 


Since M, = 1, the mean number of pulses per beam position in the noise- 
only situation is N= 1 + a,Np for a system with a single-pulse stage A, and 
INS 18K) 2 a,aNp for a system with a sequential stage A. 

For each Np there is a value of N for which Pp is maximized when the SNR 
and FAR are fixed and optimum thresholds are used. The value of N at which this 
maximum Pp occurs is close to 1 for small Np but in typical cases reaches 
N= 1+Np for Np greater than about 6 or 8. As N is increased while holding S, 
a,ap_/N and Np constant, a, increases, and therefore fh, decreases and (since ap 
decreases for fixed a 40p/N) Bp increases. For the sequential case the increase 
of a, with N depends on the fact that a, increases with E,(n) (see fig. 7) for tests 
(fixed C) of the class considered. The effect of the opposite behaviors of the 
per-stage miss probabilities on Pp) = (1-84) (1-fpg) accounts for the existence 
of a maximum Py. Graphs of Pp versus N are given as figure 8 for a system 
with fixed-sample stages and as figure 9 for a system with a sequential first- 
stage. (Rice distribution; S= 8 dB/pulse; M,= 1; a,4ap = N x 107°; optimum 
Kp for each Np and N.) Since a larger N results in fewer opportunities to detect 
the target, the actual optimum values of N (those which maximize the cumulative 
detection probability) are smaller than those for which Pp is maximized. One 
conclusion from these results is that when Np is limited to a relatively small 
value by certain factors in a particular application — possibly by the length of 
time the target remains in the resolution cell or by high-target density, clutter, 
or other conditions occasionally present which cause a large proportion of the 
stage A’s to alarm — then increasing N beyond a certain value will not improve 
even the single-scan detection probability. In fact, N can often be set con- 
siderably smaller than that value with little loss in Pp and perhaps no loss in 
the cumulative detection probability. 


*These are obtained from equations (5:20), (5:19) and (5:23) in reference 1 by letting 


P, = 1-p, > % and logi=™ = (C-1) log = (Ga) leg == 
a 


1 Po 
of B when p(S) = p,. They are also derivable by the method of difference equations 
y9)it 


°, where B, is the value 


(‘*the classical ruin problem, or by using Markov chain theory. 


10m 


a 10-4 


10-5 


10-6 


1.0 1.1 1.2 1.3 1.4 1.5 


Figure 7. False-alarm probability vs average sample size (noise-only case) 
for sequential tests of the type employed in stage A. 
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Figure 8. Single-scan detection probability vs mean number of pulses (noise- 
only case) per beam position, for a two-stage test with fixed-sample stages. 
Rice distribution; S= 8 dB/pulse; My = 1; Nq = 1; optimum Kg for each Ng 
and N; FAR = 3.6 Mg PRF x 10°, where Mp may vary. 


Figure 9. Single-scan detection probability vs mean number of pulses (noise- 
only case) per beam position, for a two-stage test with a sequential (C= 4) 
stage A. Rice distribution; S= 8 dB/pulse; M, = 1; optimum Kg for each 
Np and N; FAR = 3.6 Mg x PRF x 10°, where Mg may vary. 
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For fixed N, FAR, and SNR, there is an optimum Np and, therefore, an 
optimum a,. The existence of an optimum Np has an explanation similar to that 
for N above: as Nz is increased while holding N, a,ap/ N, and S constant, a, 
decreases and therefore B , increases and B, decreases. Figure 10 illustrates 
this for a typical fixed-sample case (Rice distribution; S= 8 dB/pulse; Ma = 1; 
N= 1.2; a,ap = 10; optimum Kp for each Np). The optimum Np tends to 
gradually increase with N (as can be seen in figs. 8 and 9); however, because 
of the nature of binomial tests it conceivably can dip, as it did for the Nn = Z 
case in figure 6. For most sequential cases, the Np giving the maximum Pp 
is too large to be practical; while the improvement of Pp with Np is noticeable 
up to about Np = 8, further improvement is not enough to compensate for the 
disadvantages (when clutter or many targets are present, etc.) of having long 
second stages. 

The value 1.2 is used for N in many of the figures primarily because it 
is small enough that rough comparisons with the N = 1 single-stage system are 
possible and also because the higher Pp that would be achieved by using a 
larger N might not compensate for the smaller average number of scans across 
the target and for the longer second stage required for optimality. 


Figure 10. For fixed M,, N, FAR, and SNR, the manner in which the per-stage 
miss probabilities change with Np produces an optimum Np. 


Even when the number Np of stage-B pulses is not constrained by other 
considerations, it can be optimized for at most only a few adjacent range bins 
for a target of given cross section. For the particular case graphed in figure 11 


OPTIMUM Nz 


(Probability of detecting 
design target if the 


Ne optimum for that 
iB OP 
particular range 
is the one used.) 


Figure 11. Single-scan detection probability vs normalized range, for a target 
of design size. 


(fixed-sample stages; Rice distribution; S, = 0 dB/pulse; S(r) = —40 log,,r; 

MW y= Ig INC = 1528 %&4 Op, = 107°; optimum Kp for each Np), the dashed curve is 
a bound for any curve with Np fixed — a curve for an Np somewhat larger than 12 
would be closer to the dashed curve only at outer ranges. Probably the more 
practical procedure is to optimize the test for an SNR that corresponds to a 
chosen value of Pp, rather than to optimize it for the SNR that corresponds to an 
arbitrarily chosen range and design cross section. Figure 12 (fixed-sample 
stages; My = ie = 12 optimum Kp for each point; Pp varied by changing the 
SNR) shows how the optimum Np changes considerably with Pp: The correspon- 
ding optimum values of the stage-A false-alarm probability (optimum a iy = a2 
[optimum Nz when N = 1.2) for these cases are plotted in figure 13. These 


27 


Optimum Np 
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Figure 12. Optimum number of stage-B pulses vs single-scan detection probability. 
Optimum Kp for each point; M, = 1; N= 1.2; Pp varied by changing the SNR. 
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Figure 13. Stage-A false-alarm probability (corresponding to optimum Np in 
figure 12) vs single-scan detection probability. 


results indicate that the percentage of times a second stage is necessary in the 
noise-only case is between about 0.5 percent and 4 percent when the two-stage 
test is optimized for a Pp between about 0.2 and 0.8. It can be seen in figures 
14 and 15 that the ratio of per-stage miss probabilities of an optimized test tends 
to fall between about 0.04 and 0.1 (a large variance results from the binomial 


oe) 


Bp/Ba 


0.001 


30 


RAYLEIGH 


0 0.2 0.4 0.6 0.8 1.0 
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Figure 14. Ratio of per-stage miss probabilities vs single-scan detection 
probability for the optimized tests of figures 12 and 13. 


and small-sample nature of the tests) even as the SNR (and therefore Pp), the 
FAR, the underlying distribution, and N vary. (Fig. 15 is for fixed-sample stages, 
Wig = lp INA Selly iy Che = IN x 10°°, optimum Np and Kp per point. Successive 
segments correspond to increasing values of Np; the ratio increases fairly 
smoothly with N for constant Np.) 


—O— RICE; S= 8dB/pulse 
——A%e- RAYLEIGH; S= 11 dB/pulse 


1.0 1.2 1.4 1.6 1.8 2.0 
N 


Figure 15. Ratio of per-stage miss probabilities vs mean number of pulses 
(noise-only case) per beam position, for optimized tests. Optimum Np and 
Kp for each point; My = 1; Ny = 1; «4 ag = N x 10°. 
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Comparison of Systems 


In figures 16 through 19, single-scan detection probabilities are plotted 
for a single-stage system, a two-stage system with fixed-sample stages, and two 
cases of two-stage systems with sequential first stages. In each figure, the two 
sequential cases (C= 4, Np = 10) and (C= 6, Np = 8), require approximately 
the same average dwell-time in the signal case, and for large signals the values 
of N(S) for the sequential cases and for the fixed-sample case are all approxi- 
mately 13. A comparison of figure 16 with figure 17 will show that the power 
savings of the two-stage systems over the single-stage system increase if the 
FAR is reduced by decreasing a, %p anda. Figures 16 through 18 do not allow 
a completely fair comparison of the single-stage case with the two-stage cases 
because of the larger number of opportunities to detect the target in an N= 1 
single-stage case than in N = 1.2 cases, so a comparison for N = N= 2 was 
included as figure 19. 

As was seen in figure 7, the sequential test with C = 6 is capable of 
providing a smaller a, than the C= 4 test with approximately the same E,(n). 
However, since a very small a, isnot efficiently utilized in a two-stage system, 
the detection probability of a system employing the C = 6 test is only slightly 
better than that of a system employing the C = 4 test and the same Np. The 
C= 4 test is generally a good choice since tests with larger C result in larger 
sample sizes in clutter or signal situations. The tendency of the sequential 
system toward a long first stage when clutter is present or the signal is small is 
evident in figure 20. 
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Figure 20. The mean number of pulses in a sequential stage-A vs the probability 
that the receiver output voltage exceeds the quantization level. 
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MULTIPLE-RANGE-BIN FIRST-STAGE 


In most of the practical applications of two-stage systems, it will be 
necessary to use at least several first-stage bins to cover the search zone. 
Figure 21 presents plots of Po vs S for a Rice distribution for various combinations 
of M, and N, where FAR = 3xPREFxM,Mp,~x10™. As mentioned in the dis- 
cussion of figure 3, for a fair comparison of cases having different values of 
M,, the M, bins must cover the same search zone and M, Mp must have the 
same value. From a comparison of these figures the following observations 
can be made. 


1. An increase of about 0.6 to 1.0 dB/pulse in the SNR (the SNR required 
for a single-scan detection probability of 0.5) is necessary as M4 is quadrupled 
(under the assumptions made in the earlier discussion under ‘‘Resolution’’). This 
same observation was mentioned in connection with figure 3. The amount of this 
increase tends to be smaller for the larger values of M,. 


2. For fixed N (and Pp = 0.5), the decrease in the optimum value of Np 
with an increase in M, is often considerable. 


3. For fixed M,, the optimum value of Np increases with N; this is con- 
sistent with the results for M, = 1 shown in figure 6. 


4. The optimum value of Np decreases with an increase in Pp (equivalently, 
with an increase in S), as was shown for the M, = 1 case in figure 12. 


One important conclusion drawn from these figures is that often an Np 
considerably smaller than the optimum Np can be used without significant loss 
in detectability; this is fortunate, since a small value of Np is desirable in 
situations involving occasional clutter or high target density. 


S (in dB/pulse) 


Figure 21. Single-scan detection probability vs SNR, for various 

oe —5 
values of N and M,. Rice distribution; FAR = 3 PRF x Ma, Mp * 10”, 
where Mp, may vary. 
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Figure 21 (Continued) 
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Figure 21 (Continued) 
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Figure 21 (Continued) 
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Figure 21 (Continued) 


S (in dB/pulse) 
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Figure 21 ( Continued) 


S (in dB/pulse) 


APPLICATIONS INVOLVING CHANGES IN ENERGY 


When employing a two-stage procedure with any particular radar, the full 
power of the radar must, of course, be used in both stages if the best possible 
detection probability is to be achieved. This policy generally results in using 
the same per-pulse energy in both stages if the radar is average-power limited 
but is not (for the waveforms used) peak-power limited. If continual operation 
at maximum power is impractical, as when the source of energy is so limited by 
its cost or weight that only a certain average energy can be transmitted per beam 
position, then the optimum division of energy between the two stages (equivalently, 
the optimum per-pulse energy for each stage) must be determined for each detec- 
tion situation under consideration. The determination of these optimum energy 
levels may involve boundary conditions imposed by peak-power and/or average- 
power limitations. 


CONCLUDING REMARKS 


Two-stage procedures promise considerable saving over conventional 
procedures when the target density is low and the radar video is clutter-free. In 
other environments, two-stage detection systems, like all sequential detection 
systems, spend too much time in too many sectors. This can be alleviated to 
some extent by automatically switching to a fixed dwell-time mode in cluttered 
or jammed sectors. 

For cases where the average number of pulses per beam position is less 
than about two and a single range-bin is examined in the first stage, calculations 
showed that the detector with a sequential first stage has a slightly higher de- 
tection probability than the detector with fixed-sample stages; however, this 
saving hardly seems enough to compensate for the more complex implementation 
required by the sequential system and for its slowness in clutter situations. For 
applications which call for a large average number of pulses per beam position, 
it seems likely that a detector employing a carefully designed sequential test in 
the first stage would have a substantially better detection capability than a de- 
tector with fixed-sample stages, especially if only one or a very few range bins 
are used in the first stage. The sequential test becomes less suitable as the 
number of first-stage range bins increases since, as a number of people have 
demonstrated for a variety of multiple-resolution-element sequential tests, the 
saving of sequential tests over fixed-sample-size tests decreases as the number 
of resolution elements increase. This latter property of sequential tests is the 
reason for considering their incorporation into two-stage detection procedures 
when high resolution is required; the benefits of sequential detection may be 
realized in the detection function of the coarse-resolution first stage, while the 
fine-resolution requirements are fulfilled in the second stage. 


45 


REFERENCES 


1. Wald, A., Sequential Analysis, Wiley, 1947 


2. Westinghouse Research Laboratories Report 412FF512-R2, Analysis of a Two- 
Stage Signal Detection System, by C. W. Helstrom, 22 April 1960 


3. Finn, H. M., ‘‘A New Approach to Sequential Detection in Phased Array Radar 
Systems,”’ p. 4-3 — 4-12 in National Winter Convention on Military 
Electronics, Conference Proceedings. 2d, Professional Group on Military 
Electronics, Institute of Electrical and Electronics Engineers, Los Angeles 
Section, 1963 


4. Rome Air Development Center Report TDR-64-122, Optimum Radar Surveillance 
Mode Study, by H. M. Finn, July 1964 


5. Finn, H. M., ‘‘Energy-Variant Sequential Detectors for Fluctuating Targets,” 
p. 444-468 in Rome Air Development Center Report TDR-64-225, v. 1, 
Proceedings of Symposium on Electronically Scanned Array Techniques and 
Applications, July 1964 


6. Rome Air Development Center Report TR-65-77, Optimum Radar Surveillance 
Mode Study Program, by H. M. Finn, September 1965 


7. Marcum, J. I. and Swerling, P., ‘‘Studies of Target Detection by Pulsed Radar,’’ 
Institute of Radio Engineers. Transactions: Information Theory, v. IT-6, 
p. 59-308, April 1960 


8. Rand Corporation Memorandum RM-4338-PR, A Two-Step Sequential Procedure 
for Improving the Cumulative Probability of Detection in Radars, by 
L. E. Brennan and F. S. Hill, Jr., October 1964 


9. Brennan, L. E. and F. S. Hill, Jr., ‘A Two-Step Sequential Procedure for 
Improving the Cumulative Probability of Detection in Radars,”’ Institute 
of Electrical and Electronics Engineers. Transactions: Military Electronics, 
v. MIL-9, p. 278-287, July/October 1965 


10. Baum, R. F., ‘‘Simple Sequential Detection Tests,’’ Institute of Electrical and 
Electronics Engineers. Proceedings, v. 52, p. 1251-1252, October 1964 


11. Feller, W., An Introduction to Probability Theory and Its Applications, v. 1, 
2d ed., Wiley, 1957 


UNCLASSIFIED 


Security Classification 


DOCUMENT CONTROL DATA-R&D 
(Security classification of title, body of abstract and indexing annotation must be entered when the overall report is classified) 
1. ORIGINATING ACTIVITY (Corporate author) 2a. REPORT SECURITY CLASSIFICATION 


Naval Electronics Laboratory Center UNCLASSIFIED 
San Diego, California 92152 2b. GROUP 


3. REPORT TITLE 


TWO-STAGE AUTOMATIC RADAR DETECTION 


4. DESCRIPTIVE NOTES (Type of report and inclusive dates) 
Research Report, January 1965-August 1967 


5. AUTHOR(S) (First name, middle initial, last name) 


R,. A. Worley 


6. REPORT DATE 7a. TOTAL NO. OF PAGES 7b. NO. OF REFS 
7 May 1968 48 11 


8a. CONTRACT OR GRANT NO 9a. ORIGINATOR'S REPORT NUMBER(S) 


b. Prosectno. SH O01 02 05, Task 6072 


3 
(NELC D11771) ips 


9b. OTHER REPORT NO(S) (Any other numbers that may be assigned 
this report) 


c. 


10, DISTRIBUTION STATEMENT This document is subject to special export controls and each 
transmittal to foreign governments or foreign nationals may be made only with prior 
approval of Naval Electronics Laboratory Center, San Diego, California 92152. 


11. SUPPLEMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 


Naval Ship Systems Command 
Navy Department 


~ oie ea Multiple-stage detection systems which allow stage-to-stage changes in 
the radar characteristics and in the statistical decision process are described, and 
the literature on such systems is briefly reviewed. A binomial two-stage detection 
system with coarse-fine resolution is then outlined and the optimization of some of 
its design parameters considered, with appropriate curves given. Some special 
cases of these binomial two-stage detectors are investigated in detail; for these, 
detection probabilities and optimum design parameters are plotted as functions of 
other design parameters and of the signal strength, and comparisons made with a 
conventional detection method. 


DD Ve ..1473) Ace UNCLASSIFIED 


S/N 0101-807-6801 Security Classification 


UNCLASSIFIED 


Security Classification 


KEY WORDS 


Radar targets - Detection - 
Mathematical analysis 


DD ..1473 (sack) UNCLASSIFIED 
(PAGE 2) Security Classification 


CHIEF OF NAVAL MATERIAL 
DNL 5 
COMMANDER, 
SHIPS 035 
SHIPS 2052 
SHIPS 204113 
COMMANDER, NAVAL AIR SYSTEMS COMMAND 
AIR 5104 
AIR 5330 
AIR 604 
COMMANDER, 
ORD 0322 
ORD 083 
ORD 9132 
COMMANDER, 
ELEX 913 
COMMANDER, 
CODE 6170 
CODE 6175.02 
CODE 61758 
CODE 6176 
CODE 6179B 
CHIEF OF NAVAL PERSONNEL 
PERS 118 
CHIEF OF NAVAL OPERATIONS 
OP-03EG 
OP-345 
OP-07T 
OP-701 
OPp-724 
OP-09B85 
op-922Y4C1 
CHIEF OF NAVAL RESEARCH 
CODE 418 
CODE 455 
CODE 461 
COMMANDER IN CHIEF 
US PACIFIC FLEET 
CODE 93 
US ATLANTIC FLEET 
COMMANDER OPERATIONAL TEST AND EVALUATION 
FORCE 
DEPUTY COMMANDER OPERATIONAL TEST AND 
EVALUATION FORCE, PACIFIC 
COMMANDER CRUISER-DESTROYER FORCE 
US PACIFIC FLEET 
CODE 425 
US ATLANTIC FLEET 
COMMANDER SUBMARINE FORCE 
US PACIFIC FLEET 
COMMANDER ANTISUBMARINE WARFARE FORCE 
US PACIFIC FLEET 
US ATLANTIC FLEET 
COMMANDER FIRST FLEET 
COMMANDER SECOND FLEET 
COMMANDER TRAINING COMMAND 
US PACIFIC FLEET 
US ATLANTIC FLEET 
COMMANDER AMPHIBIOUS FORCE, 
FLEET 
COMMANDER SERVICE FORCE, US ATLANTIC FLEET 
COMMANDER, DESTROYER DEVELOPMENT GROUP, 
PACIFIC 
COMMANDER FLEET AIR WINGS, ATLANTIC FLEET 
NAVAL AIR DEVELOPMENT CENTER 
LIBRARY 
NAVAL MISSILE CENTER 
TECHNICAL LIBRARY 


NAVAL SHIP SYSTEMS COMMAND 


(2) 


NAVAL ORDNANCE SYSTEMS COMMAND 


NAVAL ELECTRONIC SYSTEMS COMMAND 


NAVAL SHIP ENGINEERING CENTER 


US PACIFIC 


CODE N3232 

PACIFIC MISSILE RANGE 
CODE 3250 

NAVAL AIR TEST CENTER 
NANEP 


WEAPONS SYSTEMS TEST DIVISION 
NAVAL WEAPONS CENTER 
CHINA LAKE 
CODE 753 
CORONA LABORATORIES 
TECHNICAL LIBRARY 
FLEET COMPUTER PROGRAMMING CENTER 
PACIFIC 
CODE 21 
ATLANTIC 
TECHNICAL LIBRARY 
NAVAL WEAPONS LABORATORY 
LIBRARY 
PEARL HARBOR NAVAL SHIPYARD 
CODE 246P 
PORTSMOUTH NAVAL SHIPYARD 
CODE 242L 
SAN FRANCISCO NAVAL SHIPYARD 
HUNTERS POINT DIVISION 
PHILADELPHIA NAVAL SHIPYARD 
CODE 249C 
CHARLESTON NAVAL SHIPYARD 
CODE 240 
NAVAL SHIP RESEARCH & DEVELOPMENT CENTER 
CARDEROCK DIVISION 
LIBRARY 
NAVY MINE DEFENSE LABORATORY 
CODE 716 


INITIAL DISTRIBUTION LIST 


NAVAL TRAINING DEVICE CENTER 
TECHNICAL LIBRARY 
CODE 02 (2) 
NAVY UNDERWATER SOUND LABORATORY 
LIBRARY 
NAVAL CIVIL ENGINEERING LABORATORY 
L54 
NAVAL RESEARCH LABORATORY 
CODE 2027 
CODE 5250 
CODE 5320 
CODE 5330 
BEACH JUMPER UNIT ONE 
BEACH JUMPER UNIT TWO 
FLEET ASW SCHOOL 
TACTICAL LIBRARY 
NAVAL UNDERWATER WEAPONS RESEARCH AND 
ENGINEERING STATION 
LIBRARY 
OFFICE OF NAVAL RESEARCH BRANCH OFFICE 
BOSTON 
CHICAGO 
LONDON 
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